Quantum interference and local field effects in a medium of V-type three-level atoms 
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We study the intrinsic optical bistability displayed by a small sample of V-type three-level atoms 
induced by the near dipole-dipole interaction. The use of the coherent state properties in the 
limit of the generalized second-order Born approximation for BBGKY-hierarchy of equations for 
the reduced density operators allows one to derive the operator describing the near dipole-dipole 
interaction (local field correction) with interference terms involving no additional assumptions. The 
dynamics of populations of the excited states and the total spontaneous intensity are analysed as 
functions of the external laser field strength allowing for the off-diagonal structure of the local field 
and relaxation operators. 
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I. INTRODUCTION 

Effects arising during resonant light scattering in dense 
optical media have been posing both fundamental and 
practical interest for modern quantum optics and its ap- 
plications. One of the most studied range of phenom- 
ena in this aspect is linked to the interaction of closely 
positioned atomic systems through their radiation field 
and combines a number of well investigated cooperative 
effects This type of phenomena includes local-field 
effects or the near dipole-dipole interaction, a striking se- 
quence of which is intrinsic optical bistability predicted 
and studied theoretically 0, EJ d 13 and observed in ex- 
periments with impurities embedded into crystalline ma- 
trices or doped glasses 0, Q ■ 

As it was demonstrated by the recent studies con- 
sideration of the near dipole-dipole interaction appears 
to be significant when investigating effects induced by 
quantum coherence and interference. The latter ones 
lead to such well known phenomena as elect romagnet- 
ically induced transparency 8), gain without inversion 
@ , quenching of spontaneous emission 0, 0] , effect of 
quantum interference on its angular distribution |l2j | and 
resonance fluorescence |l3j . 

In most primary theoretical works the bistable re- 
sponse of a system of two level atoms was studied by 
means of the generalized Maxwell-Bloch equations where 
the effective field driving an atom was introduced using 
the Lorentz-Lorenz relation [j| . Since this correction de- 
pends on the value of induced polarization of atoms, it 
gives rise to a nonlinearity (feedback condition) which 
produces multistability and specific hysteresis loops. 

Later, owing to the fact that the bistable behaviour 
was observed in the optical response from dope atoms 
placed in a dielectric host, a number of suitable theories 
were developed but all still exploited the original princi- 
pals and approaches 0,0- However, it has been demon- 
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strated that an atom in a dielectric medium can have 
a notably different spontaneous decay rate 0, 0, 0] . 
Also, for these case the strict condition forbidding real- 
ization of such bistability mechanism cased by collisional 
broadening was eliminated 0, 0] . 

The local field correction in the equations describing 
multilevel systems produces some nonlinear cross terms. 
The local field effects on electromagnetically induced 
transparency was studied in [Tjj, the same influence on 
lasing without inversion was considered in J2jj. In [2l| 
they report the effect of increasing the index of refrac- 
tion in an absorptionless medium and a mechanism of 
piezophotonic switching between the absorption and gain 
modes of the system. The local field effects in the mech- 
anisms of interference were studied using the tree-level 
model of A-configuration. The cascade type was analyzed 
in [221 ]. It is worth noting that quantum interference 
has been also considered as an alternative mechanism to 
reach bistable behaviour in three-level systems in cavi- 
ties [23| . The interest to the V-type model was displayed 
while studying energy level splitting with quantum co- 
herence and local field ,24|. 

Turning back to the issue of intrinsic bistability it must 
be noted that some recent experimental results leave the 
area open to searching for new interpretations of its na- 
ture [25| and suggest alternative approaches for descrip- 
tion of observed phenomena |2y, |23 • 

The aim of this work is to create a fully quantum- 
mechanical approach to allow for the local field correction 
in the laser field interaction with a dense medium of mul- 
tilevel atoms. Based on the hierarchy of equations devel- 
oped by Bogolubov, Born, Green, Kirkwood, Yvon, et a!., 
(BBGKY - hierarchy) for the reduced density operators 
treated in the generalized second Born approximation we 
will get a self-consistent set of atomic and field equations 
and equations for their correlations. With the use of co- 
herent state theory we will derive the terms describing 
the near dipole-dipole interaction containing off-diagonal 
interference terms providing specific bistable behaviour. 

This work is organized as follows. In Sec. [H]we derive 
the self-consistent set of equations including the equa- 
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tions describing evolution of atomic and field subsystems 
and their correlations. We get the explicit expression for 
spectral and total spontaneous intensity. In Sec. If HI the 
results of our numerical analysis are reported. The dy- 
namics of populations of the excited atomic states, as 
well as atomic polarization characteristics and total in- 
tensity of resonance fluorescence is studied. Finally, in 
conclusion we summarize the results obtained in our re- 
search. 



II. HAMILTONIAN AND BASIC EQUATIONS 

We consider a collection of N three-level atoms inter- 
acting with a quantized electromagnetic field in a vol- 
ume of linear dimensions much smaller than the resonant 
wavelength. The atomic level configuration is chosen to 
have closely separated upper excited level |2) and lower 
excited level |3) both coupled to a common ground state 
|1). The Hamiltonian which describes the system in the 
dipole approximation can be written in the frame rotat- 
ing with the laser frequency ujl and in units of 1/h as 



/v 



JV 



k k 

3=2,3 

V AF = J2 iX * (S+a k -alST), 



(1) 



where and ak are the creation and annihilation opera- 
tors for the field mode with wave vector k and frequency 
o>k, detuned from the laser field frequency: v k = Wk-wj,. 
Similarly, Aj — ujji — ujl is the detuning between the fre- 
quency LUji of the atomic transition \j) — > |f) (j = 2,3) 
and the laser frequency. Sjj — = and 

Sj~ = are the atomic projection operators for atom 

j. The coupling constant between the field mode and the 



atomic dipole transition = \\pji 



where 



Mil 



k ;'^"v nw 

is the matrix element of the dipole transition and W is 
the quantization volume. In the small-volume limit we 
neglect the spatial distributions and, thus, consider the 
system's parameters to be scalars. 

The interaction of the field with the collection of atoms 
we treat by means of the atom-field density operator p 
which obeys the equation 



i~P = [2M> 

JV 



P — Patom <K> P field 



k 



(2) 



Here pa and pf denote a single atom and a single field 
mode density operators respectively. 

Let us assume that initially the system consists of 
two non-interacting subsystems that are the collection 



of atoms and the laser field. The laser field density op- 
erator one can represent in terms of coherent states |«k) 
psl ] and write the corresponding initial condition as: 



JV 



p(o)=n^(°)ni a k>Ki 



(3) 



Further, in order to obtain self-consistent atomic and 
field equations we will follow the approach developed in 
papers The first necessary step it suggests 

is to eliminate the laser field from the above condition 
© which will then yield the density operator describing 
evolution of the atomic system and the scattered electro- 
magnetic field. To accomplish this we use the unitary 
transformation: 



p coh x 



p — exp(iffpi)(J^ D(ak)) exp(— iHpt) 

k 

x exp(iffirt)(]^J Z? + (ak)) cxp(— iHpt), 

k 

where -D(ak) is a coherent state operator |28j : 
D(a k ) = exp(a k a k ~ a t a k.), -D(«k)|0k) = l«k) 



(4) 



(5) 



The proper Hamiltonian in Eq. for p coh can be writ- 
ten explicitly using the well known properties of D{a\ 5 ) 
M: 



D + (a\ c )al c D(a\ c ) = + a£, 



(6) 



and expressions 

exp(zu;kak a k^) a k exp(— iujua^akt) = ak cxp(— iwki), 
exp(ia>ka k aki)a k ex P(~ i^ k a k aki) = a ik ex P(* w ki), 



in the form: 

H coh = cxp(-iHpt) 



(7) 



(II ^ + K)): 



x exp(iH F t)(H — Hp) exp(-iHpt) x 



:(Y[D(a k ))+H F 



k 
JV 



exp(iHpt) 



(8) 



JV 



E^+E^+EE 1 ^ 



where 



vtf) 



V* = E E ^ i^ s i 

k .7=2,3 

= E^ 

.7=2,3 



-V A 



) = 



(9) 



One can see that Va indeed appears to be the conven- 
tional term which describes the semi-classic interaction of 
three-level atoms with an external laser field where V,- is 
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the corresponding Rabi frequency while the other terms 
remain the same as in Yet, the initial condition in 
this picture transforms into 



where 



N 



5,i (o)=n< i (°)ni°i<>< o ki, 



(10) 



where the field part is in the vacuum state. Consequently, 
by solving Eq. ((2J) in this representation we can get the 
field density operator which contains information about 
the scattered field only. 

Our following treatment of the system is based on 
the BBGKY-hierarchy of equations for the reduced den- 
sity operators [29L |32(. Let us remind that we con- 
sider the system which contains two sorts of particles, 
i.e., atoms and radiation field photons, so the notations 
below are indexed as A and F identifying either their 
atomic or field nature. Primed indices refer to other par- 
ticles of the same sort. Moreover, we assume that no 
interaction between particles of one sort occurs, that is 
Vaa' =Vff' = 0. For the reduced density operators we 
drop the upper index coh as we consider them everywhere 
in representation Q until otherwise stated. The equa- 
tions of the BBGKY-hierarchy one can write in terms of 
two and three particle correlation operators gAF, Qaff 1 
and gAFA' defined as [13 



PAF = PAPF + gAF, 
RAFF' = PAPFPF' + gAFPF' 
PAF A 1 = PAPF PA' + gAF PA' 



' 9AFF' i (11) 
gAFA' , 



where [...] denote the terms obtained from permutations 
of particle indices in the product of correlation and den- 
sity operators. The first four such BBGKY-hierarchy 
equations of our interest are 

d 

% ~dt PA ~ \- Ha,PA ^ = Tr F'[yAF',gAF'] 1 



^~alPF - [H F ,p F ] = Tr A ,{V A > F ,g A ,p], 

d 9 

igj-dAF - [H AF ,g AF ] = [V AF ,p A p F ]+ 
+L AF + IL AF + C A > F >, 
L AF = [V AF , g AF ], 

U AF = Tr A ,[V A > F ,g AA ,p F ] +Tr F '[V AF ',gFF'pA], 
Ca'F' =TrA'[VA'F,gAFA'] +Tr F '[V AF >,g AFF >], 
d 

i-^9AFF' - [H AFF ,,g AFF '\ = 

= Iaff' + L AFF i + Haff' + C A " F ", 
Iaff' — [Vaf + V AF ' , p A p F p F >], 
L AFF ' — [V AF + V AF ' , g AFF >], 
Uaff> - 0(4), 
C A ,, F ,, = 0(4). 

Here we use the effective Hamiltonians: 



(12) 



H A = H A + U A , 
H F = H F + U F , 

H F + U A + U F , 
H A + Hp + Hp, + U A + U F + Up, 



Haf " Ha 



(13) 



H° 

n AFF 



U F = Tr A V AF p Al 
U A = Tr F V AF pp. 



(14) 



This set of equations contains two particle correlation op- 
erators g AA ' and <?ff', however, the approximation to be 
used ignores their contribution. Thus, it can be already 
regarded as a self-contained system. Namely, we treat the 
hierarchy in the generalized second Born approximation 
[32| which implies 



Iaff' — > [V AF >, p F ,g AF ], 
L AF = U AF = L AFF , = U AFF , 



0. 



(15) 



Besides, we consider contribution of the only three par- 
ticle correlation g AF F , . Finally, the equations describing 
the system read as follows: 

d 

l ^lP A ~ \ h a, Pa] =Trp,{V AF ,,g AF ,}, 
f 

i~QlPF ~ [Hf,Pf] = Tr A ,{V A ,p,g A >p}, 



igidAF ~ [H AF ,g AF ] = [V AF ,p A p F ]- 



(16) 



l di 9AFF ' 



+Tr F ,[V A p,,g AF p>], 
[H AF pi,g AF p'] = [V AF >, p F ig AF \. 



The initial conditions are 



p A (0) = p% p F (0) = |0 k )(0 k |, 
5af(0) = 0, g A FF'(0) = 0. 



(17) 



In this model we can consider the radiation to be leaving 
the interaction volume irreversibly or, in other words, we 
can ignore multiple scattering of spontaneously emitted 
photons. Thus, in further transformations of the right 
hand side terms in Eqs. (|16|) we use the first approxima- 
tion for the field density operator and treat it as 



PF = |0 k )(0 k |. 



(18) 



It is now advantageous to carry out the analysis in the 
wave picture for which one needs to rewrite Eqs. ll(j|) 
substituting 



p = e -iH F t~ e iH F t 



iH F t 

+ H.c. 



V AF = c z V AF 6~ 

.7=2,3 _ 

H A = H A + U A , 
Hp = Up, 
H° AF = H A + U A + U F , 
H AF p,=H A + U A + Up + Up,. 

The Up potential is explicitly 

JV 3 

U F = Tr A V AF p A = } j } j (n\V AF p A \n) 



(19) 



N 

3=2,3 



(20) 



X^mpA^e-^+H, 
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For the volume smaller than the wavelength it is appro- 
priate to take the sum over atoms to multiplication by the 
number of the particles N. By introducing the notation 



t k (t) = -Nj2^ ) (j\p\l)e iVkt , 
we reduce (7f to 

uf = *(Ck(*)4 - ciWok), 



(21) 



(22) 



where U F potential is, indeed, the quantum representa- 
tion of field induced polarization of the matter [29] . One 
can see that the operator exponent exp(— iU F ) is by def- 
inition a coherent state operator [2^ and, therefore, the 
sequent relations 

T{t) = exp(-iU F )(t) = exp(-ty(t))£>(&(t)), 
^/3k = Ck, ^0=-z/2(/3 k C k -/3 k Ck), 

are valid. Now we transform to a new representation by 
defining 



p = T(t)p p T-\t) 



(24) 



The interaction operator in the right hand terms of Eqs. 
()16J) is modified according to the properties of the coher- 
ent state operator similar to JJJJ: 



e X p(^(t)) J D+( ( 5 k (t))a k exp(-i0(i))£)(/3 k (t)) 

= (a k + &(«)) 
exp(^(t))D+(/? k (t))4 exp(-i4>{t))D(p k (t)) 

and explicitly takes the form 

T- 1 VafT = Vaf + V£ f , 
Vl F = £ (/3 k 5+ e"^' - /3 k 5 7 - e 

i=2,3 



(25) 



(26) 



Correspondingly, the U A potential in the Hamiltonian 
transforms into 

T- X U A T = U A + U£, 
Ul = Tr F V£ F p F = EE< nk I^P?K) 

k «k I07\ 

U) (R. Q+ „~iv k t_ \ Z <> 



k n k 



while the effective Hamiltonians in this picture are 

El = #V = H° AFF , = H A + U A + U%. (28) 



Since contains no field operators taking the trace in 
(|27|l by definition eliminates the density operator. Now 
the expression for U? potential can be significantly sim- 
plified. Changing the sum over k into an integral and 



substituting the formal solution for /3 k (i) obtained from 
Eqs. (23J! and (|3TJ to J2ZJ we get 

U% = -iN [ S+UW) I r+ j dt>dv k + H.c, 

j,m=2,3 J ° J "k 

(29) 



where 



A j(i) A (>") e M*-*') 1 
(30) 



Next we apply the Wigner-Weisskopf approximation 
and find that the operator U p is 



i/2iV J2 



S+(j\p\l)+H.c. 



j } m —2,3 



(31) 



where = 7^ and r mj - = p^/^WTj with 7^ being the 
spontaneous emission rate from the jth excited state 
given by jj — 40/^ \ \ 2 /3hc 3 . Parameter p arises from 
integrating the non-diagonal terms (|30l) and reflects the 
fact that one needs to account for the dipole vector ori- 
entation p = P21 ■ M3i/||M2i||||M3i||- We will consider the 
cases when p takes the values 0,1, and —1 depending on 
wether the vectors are orthogonal, parallel or antipar- 
allel ^(J. It must be noted that the strict and correct 
introduction of the parameter p requires consideration of 
the system subject to degeneracy of magnetic momentum 
projection |12lll3| . 

In its turn in the limit of the approximation (|18fl the 
impact of Ua vanishes 



Ua = 2222(nu\VAFp F \nk) 



EE' 



nk\a k p F 



aLp F \n^) 



0. 



(32) 



Transforming back from the wave picture we finally get 
. d 







l ~^P A ~ \ h a, Pa] 

■ P rfj Pi 
l dt Pp ~ [ F,Pp \ 



Tr F > [V AF > 
: Tr A > [Vaf 



' Vaf'iQaf'. 
Vaf, 9a'f], 



i di 9AF ~ 9 ^ = ^ Vaf + V af^PaPf] 

+Tr F , E k [VAF' + V£ F „gP FF ,], 
dp - P ~ P P P 

i dt gAFF ' ~ [ H aff>,9aff>] = [ Vaf " + V af',Pf>9af], 

(33) 



H F 

U AF 
n AFF' 



= H A - 

= Hp, 

= H A - 
= H A 



Hp - 



U A > 

-H F , 



(34) 



U F . 



The term in the right hand side of the first equation 
describes relaxation of the atomic system. In order to 
obtain the explicit form of this operator the formal solu- 
tion for g AF (t) must be substituted. It works out if we 
temporarily change to a picture where 



9af 



PA = C(t)p A C~ 1 (t), 



- iH ^T(t) g p F T-\t) e lH "\ 



(35) 
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with C(t) and T(t) being the fundamental solutions 

d d 
for equations i^PA - [H A ,p A ] = and *^Saf ~ 

[H A ,g AF ] = respectively. 

The first and the third equations in l|33[) with no ac- 
count for contribution from the three particle correlation 
transform to 

i^p A = Tr F ,[C-\t)(V AF , +v£ F ,)C(t),C- 1 (t)x 
xT{t)g p AF ,T(t)^C{t)], 

^9af = [T{t)-\V AF , +V£ F ,)T(t),T(t)- 1 x 
xC{t)pC- l {t)T{t) PF ]. 

(36) 

Now, substituting the formal integral for Q AF into the 
first equation gives 

ijr,P=-i f dt' Tr F [C-\t)(VAF,(t) + V£ F ,(t))x 



xC(t),L(t) [TitT 1 {V AF ,(t') + V£ F ,(t')) x 

xT(t'),L-\t')pL{t')p F ]L- l (t)\ =T(p), [6<) 

where L(t) = C" 1 (t)T(i). In the limit of (JTSJ it is ob- 
vious that the trace operation here cancels all terms in 
the right hand side, denoted as T(p), except those pro- 
portional to akd^pF, PFOkfl^ and a^ppa^. This fact, to- 
gether with and (|26|l . straightforwardly demonstrates 
that r(p) eventually splits into two parts. The first part 
combines the terms to the second order in Ak and, as 
seen below, is the well known relaxation operator for 
the three-level atomic system. The second part com- 
bines the density-dependent terms to the fourth order 
in At and represents a correction which refers to the lo- 
cal field enhancement of the spontaneous emission rate 
[LH ITEl H^. However, this correction deserves an indi- 
vidual and thorough consideration and is not reflected in 
this paper. Thus, to the second order in the coupling 
constant is 



r(p) = 



mj 



-iW dt'ic^tt)^ £ S+T(t)T-\t')S-C(t')pL(t')L-\t)fi 

k ^° ^ j=2.3m=2,3 

+ L{t)L-\t>) £ J2 t ^pC(tT 1 S+T(t')T-\t)S-C(t) 

m=2,3j=2,3 

-C^WE £ rjrnSrTW-^t'^C-H^S+TWL-^t') 

j=2,3m=2,3 

-L{t)T-\t') £ Y, t rn J S m C(t')pL(t')T~\t)S+C(t)\, 

rn=2,3j=2,3 ' 



(38) 



r 



where we used the notations l|3()(l . Now as in l|29|l - H31(l we 

take the sum over modes k to an integral and apply the 
Wigner-Weisskopf approximation. Changing back from 
p A to p A we finally get 



T(p) = i 1/2 £ ((T ]m +T+ lJ )S m pS+ 

j: m — 2,3 



i P^m &j 



(39) 



By taking the same steps as applied to the last two equa- 
tions in 1)33(1 one can get the relaxation operator for g AF 
identical in its structure to (|39|1 . Eqs. (|33|l are explicitly: 



[H A , pa] = t (pa), 



d 
Of 

i^rPF - [ h f,Pf] = Tr A ,{V AF ,g%, F }, 
i dt 9 ^ F ~ [ H af>9af] = [Vaf,PaP F ] + F(g AF ). 



(40) 



The assumption that spontaneously emitted radiation 
leaves the volume without scattering allows us to deter- 



mine its spectral distribution by calculating the proba- 
bility of emission of a single photon to mode k. This 
value is proportional to the matrix element (lk|pjr|lk)- 
The equation to find this element can be obtained from 
Eqs. PHI) by using the approximation (JTHJ, i.e. 



^(ik|p£|i k ) = 
= -iJV £ ^{(Ok, j\g^ F \i k , i) + (i k , i|fl^|0k,i)}. 

3=2,3 

(41) 

Let us denote the following matrix elements as 

g { - ] = <0k|s^|lk), 9 (+) = (lk| 5 ^|0k). (42) 

If the terms in the last equation in i|40[) are multiplied 
from the left and the right by certain corresponding vec- 
tors we get the equations for g(~^ and g( + > and, in such a 
way, transform the system to an equivalent form suitable 
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for finding (l k |/j£|l k ) 
d 

l -Q t PA ~ [Ha, Pa] = T(pA), 

4<0k|p£|0 k ) = zN £ X^{{j\g^\l) + (l\gM\j)}, 

i=2,3 

4(lk|^|l k ) = -iiV £ A«{(j|. 9 (-)|1> + (l|<7 (+) |j», 

3=2,3 

. a 

'at' 



3=2,3 



• U (+) 

V = 



= [H A ,g 



3=2,3 



p a (o) = p% (o k | P £|o k )(o) = 1, (i k | P £|i k )(o) = o, 

ff (-)(0) = g (+)( ) = 0. 

Here we have added the equation for the matrix element 
(0 k |p^|0 k ) which we need to find the density operator 
describing the observed radiation. Integrating the field 
equations formally gives 

(0 k |p£|0 k ) = 1 - N V a£° f dtFfa, t), 

3=2,3 J t ° 

(l k |p?|lk> = AT V A^ f dt'F(v k ,t), ,44) 

3=2 3 ^° 

^K,t) = (i|5«|i) + (ib (+) li)- 

Since we ignore the field states with more than one pho- 
ton in mode k, transforming back from l|24ll and |@J| to 
the initial picture and dropping the terms to the forth or- 
der in A k the field density operator to describe radiation 
observed in direction other than that of the laser beam 
(one puts a k = 0) takes the form 



p F =D(f3' k (t))\0 k ){0 k \D+(f' k (t)) 



3=2,3 



N J2 x k } I )iik><i» 

4—0 1 JO 



(45) 



and represents nothing but the density operator of the 
scattered field. The spectral intensity of scattered radia- 
tion is defined as [2^ 

J k = Kk ?iw k (|/3 k (i)| 2 

■i-o i Jo 



-K k t ) 



3=2,3 



Here the factor and the exponential index « k is a con- 
stant which determines the mode loss on the quantiza- 
tion boundary or during signal detection or any other 
irrevocable loss of radiation. The need and mechanism 
of accounting for the loss of photons through a detector 
has been minutely discussed in paper [29j- There, they 



demonstrate that the corresponding additional term in 
the Hamiltonian, K[p] = — i ^ k (Kk/2)(a k a k /9— 2a k pa£.+ 
pa^aic), makes some new terms appear in the equations 
for the field density operator and for the atom-filed cor- 
relation operator. However, having this term in the ini- 
tial picture does not effect the steps taken to obtain the 
results of this work. Therefore, we straightforwardly sub- 
stitute this damping corrections to the final result. For 
the stationary case (i — > oo) we have 



I k =N huj k { N X 

j,m— 2,3 



(jK(«0 

k A k ,.2 



'i'k 



• X 



■ £a£>F(i*,oo)) (4?) 

3=2,3 



The first term in the sum here describes the elastic 
or Rayleigh scattering. The second term represents 
the spectral intensity of resonance fluorescence. Since 
F(fk,t') is a function of Xjp , substituting its formal so- 
lution to H47J1 would show that both spectral components 
are proportional to A^ A k • Thus, there rises an inde- 
terminacy in the model under consideration associated 
with the fact that in real physical systems the coupling 
constant is the function of the angle between the dipole 
transition vector and the wave vector k. Consequently, 
it makes sense to study a particular physical system to 
calculate the real spectrum for and consider angular and 
polarization distribution of the scattered field. Note that 
this issue we studied in paper |l3|. As for this work 
we shall not go beyond studying the integrated spectral 
characteristics of radiation. Integrating over k in lj46|l 
accurate to o(A2s/uj) yields the expression to determine 
the total intensity of scattered radiation. 

hotai oc N hu> (N ^2 T m j pjip* nl + T mj Pj m ), 

j,m=2,3 3=2,3 

(48) 

where u) = (u>2i + u>3i)/2. The Rayleigh component de- 
serves a separate comment. For simplicity we initially 
dropped the spatial dependencies. However, the Rayleigh 
scattering is a coherent process so the angular distri- 
bution of in this way scattered light in homogeneous 
medium should be close (same if elastic) to that of the 
source beam 0. Since, as it was stipulated earlier, we 
"observe" scattering from a different direction and the 
atomic system is chosen to be homogeneous and isotropic 
we get to see nothing but the resonance fluorescence de- 
termined by the last term in (14811 . 



III. NUMERICAL ANALYSIS AND 
BISTABILITY 

In this section we report the results of our numerical 
calculations based on Eqs. I03J)-(0BJ) to study the steady 
state behaviour of the total fluorescence intensity as a 
function of the resonant field strength. Some other rel- 
evant dependencies are shown. Let us remind that the 



7 



primary aim of this work is to present a new approach to 
studying resonant light scattering in a medium of multi- 
level atoms in a range of densities. Therefore, we have 
limited this report to just a qualitative illustration of cer- 
tain observable effects expected in such systems. In 
they studied the effect of quantum interference in the 
optical response from a three-level V-type atom. Based 
on the specific parameters selected for that study we are 
going to consider modifications likely to occur due to in- 
creasing density of the atomic medium. 

In order to find the total intensity of scattered light one 
first needs to know the components of the atomic density 
matrix which evolve according to the first equation in the 
set Using ©, (GU, d3J- and (33, and substituting 
shorter notations pj m for the matrix elements (j\pA\m), 
we get 

i~QlPjm = (Wjl — ^ml)Pjm — i/% J]] (rJPsm + Pjs^ms) 

s=2,3 

+ (V; - iN/2 T% P,l)Plm- 
9=2,3 

qmPlq), 

q=2,3 

d 

s=2,3 

+(V;-iN/2j2rj q P< 1 i)Pn- 

9=2,3 

- E Pjr(K*-iN/2 £ T+p 9 i), 

r=2,3 9=2,3 

Pij=P*j> 3,m = 2,3, pu + p 2 2 + P33 = 1 

(49) 

In its turn, from 148(1 it follows that the total fluorescence 
intensity is determined by the expression: 

Ispont oc N Tilo (72P22 + 73P33 + 2r 32 Rep 2 z) (50) 

In fact, disregarding the density-dependent corrections 
in the atom-field interaction terms one can see that Eqs. 
(|49|l are well known and have been used in a great num- 
ber of works including Jf!]. In this case, however, in 
addition to the cross terms describing spontaneous decay 
and giving rise to quantum interference effects there are 
similar cross terms in the the field correction. Besides, 
it is clear that the corrections are of the same nature as 
those obtained for the Bloch equations by replacing the 
local driving field, which produces Vj, with the macro- 
scopic field when the two fields are related by the Lorentz 
local-field condition Q . The new properties of the sys- 
tem follow from the analysis of the steady state solutions 
for Eqs. (|49|) which we have found numerically using 
the Newton's method. The stable solutions have been 
verified by examining the fluctuations of corresponding 
magnitudes. The total fluorescence intensity has been 
calculated with l(5U|) as a function of the driving field 
strength equivalently expressed through V = Vyi- All 
rates used herein as parameters we present in units of 72 
being the radiative relaxation rate from the upper excited 
state. The total fluorescence intensity I sp0 nt is as well 
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FIG. 1: Total fluorescence intensity (upper picture) and 
populations of excited states and cross correlation contribu- 
tions for p — 1 (lower picture) as functions of external field 
strength with normalized parameters: 73 = O.572, A2 = IO72, 
A3 = — 6O72, N = 300. Rabi frequenc ies for corresponding 
excited states are related by V13 = 1/73/72 V12. The solid 
and dashed curves represent stable and unstable states re- 
spectively. The thin dash-dot curves show the similar depen- 
dencies for the limit of low density, i.e., N = 1. 



normalized to the peak spontaneous intensity from the 
atoms excited to the upper state only, i.e., Iq oc Nfko"/2- 
FIG El As one can see in FIGHJand FIGEJfor a small 
number of atoms within the volume the populations of 
the excited states and, therefore, the total intensity ex- 
hibit precisely the same behaviour as the dependencies 
obtained from the equations laking the Local field cor- 
rection (thin dash-dot curves) both with p = and p = 1 
(the situation with p = — I is equivalent to the one with 
p = 1). However, starting with a certain number of atoms 
some cooperative phenomena in the region of weak ex- 
ternal field should occur. Due to non-linearity induced 
by the Local field operator the system's main feature is 
the existence of multiple solutions (multistability). In all 
plots, steady stable and unstable states are represented, 
respectively, by solid and dotted curves. In FIG^ an d 
FIG|21 we picture an example of cooperative behaviour 
(thick solid and dashed curves) in a dense sample which 
one can recognize as a hysteresis cycle. In the region of 
low values of the external field, from the origin to 
the curves demonstrate relatively weak growth. Upon 



0.6-, 




V/y 2 



FIG. 2: Total fluorescence intensity (upper picture) and pop- 
ulations of excited states (lower picture) as functions of ex- 
ternal field strength with p = and other parameters and 
notations same as in FIGQ 



reaching the threshold value V^. there occurs a transi- 
tion to the upper stable state where for V > V c +, es- 
pecially in the region of significantly stronger field, the 
nonlinear effect vanishes and the corresponding curves 
for N = 1 and N = 300 merge. When, going backwards, 
one decreases V starting from a sufficiently strong field 
and gets to the region below V c +, it is seen that the flu- 
orescence intensity (as well as the other dependencies) 
deviates smoothly from the non-cooperative path and af- 
ter V~ jumps down closing the cycle in the cooperative 
state. The value of V^. marks the upper limit for the 
Local field effects to take place and increases with the 
number of atoms. V~, in turn, appears to have a more 
complex dependence on N which we discuss later in this 
section. Meanwhile, it is essential to note certain changes 
in general behaviour of the system, like the ratio of pop- 
ulations of upper and lower excited states, in the cooper- 
ative region < V < V£. As one can see, for the small 
number of atoms the upper atomic state is much more 
populated in this region than the lower state and deter- 
mines the fluorescence intensity. The region where the 
populations become comparable (p = 0) or, alternatively, 
where the population of lower state becomes greater then 
that of the upper one (p = 1) is well beyond V£. For 
p = the greater number of atoms causes the popula- 
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FIG. 3: Real and imaginary parts of the off-diagonal density 
matrix elements as functions of external field strength with 
parameters and notations same as in FIG0 



tions of the excited states to even out everywhere in the 
cooperative region and, in case of a very dense sample 
(FIG0), makes the lower excited state more populated 
which governs the intensity curve. However, when p = 1, 
with expressively smaller values of pii and /O33 their gen- 
eral ratio remains the same, i.e., more atoms are excited 
to state 1 2) regardless of whether the number of atoms is 
increased. In both limits of p = and p = 1 the value 
of shifts to the region of stronger field with density 
being increased. For p — this dependence is demon- 
strated by the dashed-dotted curve in Further in FIGd! 
and FIG 0] one can see the field dependencies of Im(pij) 
and Re(pij) which determine the induced polarization. 
Here, in the cooperative region Im{p\i) displays a linear 
growth for both values of p parameter. Such a behaviour 
has been demonstrated for the cooperative hysteresis in 
a system of two-level atoms 35]. The value of Im{p\s) 
also displays considerable growth (decrease for p = 1) in 
the cooperative region. 

FIGEldemonstrates the behaviour of V^ rit and V~ it as 
functions of the number of atoms for p = 0. The curve 
describing V~ r (solid line) has ruptures occurring for the 
values of N\ and N2. In the region of Ni < N < N2 the 
function of V~ is characteristic in its large values (see 
the Fig.). It means that at a certain number of atoms 
N > Ni the hysteresis region (V c + — V~) is narrowed 
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FIG. 4: Real and imaginary parts of the off-diagonal density 
matrix elements as functions of external field strength with 
parameters and notations same as in FIG [5] 
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FIG. 5: Critical values of the field strength as functions of 
the number of atoms. 
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FIG. 6: Same as in FIG|5| (upper picture) for undercritical 
density (upper) N = 657 and critical density (lower) N = 658. 



ber of atoms the region of the stable state V < no 
longer exists. Thus, if the system is tuned so that the to- 
tal intensity of spontaneous emission is in its maximum, 
i.e., Ispont/Io is within the region of V~ < V < , one 
can increase the density of atoms slightly which makes 
the system jump (" switch" ) down to the state where the 
spontaneous emission is nearly absent. In fact the critical 
field strength transforms as V~ r — > Vj^ . A phenomenon 
similar to the one being described has been predicted 
using the model of A type three level atoms [21j. It con- 
sists in density controlled switching between absorption 
and amplification modes and is known as piezophotonic 
switching. For the wavelength of the order of 10 -5 cm 
and the value of the dipole transition matrix element 
fiji = 1CP 18 cm the switching demonstrated in FIGE1 
can occur at the density of the order of 10 18 cm" 3 and 
laser field intensity of about 1 MW-cm -2 . 



sharply and, as it is seen in the Fig., quite substantially. 
Let us see how this property effects the the observable 
quantities. FIG El demonstrates the same dependencies 
as in FIG|21but for the subcritical and postcritical val- 
ues of the number of atoms. Here, one can see that such 
a "transition" corresponds to the tangency point be- 
tween the upper stable and unstable state curves at the 
field strength having the value Vj£ . For the greater num- 



IV. CONCLUSION 

In this work we have suggested a fully quantum kinetic 
approach for studying the processes occurring during res- 
onance light scattering in dense media of V-type three- 
level atoms with account for the Local field correction. 
This approach is based on using the BBGKY-hierarchy 
for the reduced density operator. Solving the BBGKY- 



10 




V/y 2 



FIG. 7: Same as in FIG|3 (lower picture) for undercritical 
density (upper) N = 657 and critical density (lower) N = 658. 

hierarchy in the generalized second order Born approxi- 
mation we derived a self-consistent system of equations 



for the atomic and the field part of the density oper- 
ator and their correlations. The use of the properties 
of the coherent state operator in the limit of the small 
sample allowed straightforward derivation of near dipolc- 
dipole interaction (local field correction) operator in the 
system of BBGKY equations. Thus, the approach be- 
ing suggested automatically takes the proper account for 
the contribution provided by the near dipole-dipole in- 
teraction in a dense ensemble of atoms and requires no 
complementary assumptions or phenomenological proce- 
dures. Since the relaxation operator as well as the local 
field correction operator contain off-diagonal interference 
terms, the character of the intrinsic optical bistability 
displayed by the system of multilevel atoms can be sig- 
nificantly different from that of the well studied intrinsic 
bistability of two-level atoms. 

In the second section of this work we have carried out 
a numerical analysis of the equations derived and illus- 
trated the system's behaviour when both interference and 
Local field effects are taken into account (FIG^I. We 
have also demonstrated that in the absence of coherence 
between the excited states there rises a possibility for the 
density dependent switching between the steady state 
modes expressing the maximum and weak spontaneous 
intensity respectively (FIGEJ). 

We would like to thank A.N. Starostin for fruitful 
discussions and suggestions. This study was carried 
out with support of grants from the Russian Foun- 
dation for Basic Research No.02-02-17153 and No.03- 
02-06590 and grants of the President of the Russian 
Federation No.MK1565. 2003.02, No.MD338. 2003.02 and 
No.NS1257.2003.02. 
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